Spatially sequential turn-on of spontaneous emission from an atomic wave packet 
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We discuss a simple example demonstrating that spontaneous emission from "space-time- 
superposed" atomic center-of-mass wave packets is nontrivially and time-dependent modified with 
respect to the standard dipole-pattern typical of "space-superposed" wave packets. Our approach 
provides an approximate description of a nonsimultaneous interaction of electromagnetic field with 
different parts of a wave packet. 
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I. INTRODUCTION 

Standard treatments of spontaneous emission from atomic wave packets (Cohen- Tannoudji et at, 1991; Graham et 
ai, 1992; Rzazewski and Zakowicz, 1992; Ren et ai, 1992 Rzazewski et ai, 1994; Steuernagel and Paul, 1995; Stoop 
and Rzazewski, 1995) assume situations where the atomic exited state couples to the electromagnetic vacuum at all 
points of the wave packet "at t — 0" . For t > the atom starts to decay. The beginning of this process is assumed 
to take place at all points of the wave packet simultaneously. Formally this means that at t = the wave packet is in 
a superposition of the kets |r, +), where r is an eigenvalue of the atomic center-of-mass position operator, and "+" 
means an electronically excited state. What happens for t < is usually ignored. 

An experimental procedure that leads to such a situation may be based on a simultaneous excitation of the atom 
at all points occupied by the wave packet. Alternatively, one can think of situations where a metastable long-lived 
state is coupled to the vacuum reservoir by a simultaneous turn-on of an electric field in the whole region where the 
atom is supposed to reside. It is clear that the outlined scenario (below referred to as a "simultaneous spontaneous 
emission" (SSE)) does not cover all the experimental possibilities. First of all, the atomic wave packets can be 
macro scopically extended. This happens practically in all atom interferometer experiments where a phase shift is 
produced by applying different fields to different parts of the interfering wave packets. Although typical experiments 
involve static arrangements there is basically no problem with manipulating the atoms in a fully controllable time- 
depended manner. One can imagine, for example, a situation where different parts of the wave packet are either 
nonsimultaneously excited or nonsimultaneously couple to the vacuum. Actually, at the moment of writing the 
paper we are aware of an experiment where the latter possibility was practically realized. What is interesting the 
experimental data do not completely agree with the standard SSE computations. 

The works of Cohen- Tannoudji et al. (1991) and Graham et al. (1992) suggest that either there is completely no 
influence of the coherence of an atomic wave packet on the radiated photons, or that there is at least no coherence 
in the emitted light if the atomic wave packets do not overlap. The work of Rzazewski and Zakowicz (1992), on the 
other hand, shows clearly that spectral density of radiation carries some information about the shape of the wave 
packet even if one considers a two-peaked Gaussian whose peaks are separated by several wavelengths of photons and 
are essentially non-overlapping. The two-peaked wave function in position space translates into an oscillatory wave 
function in momentum space, therefore manifesting the modified probability distribution of momenta in the wave 
packet. This momentum distribution implies a distribution of Doppler shifts in the spectral density of emitted light. 
This interesting result means that, in principle, one can investigate a structure of atomic wave packets by detecting 
photons emitted during such an apparently incoherent process as spontaneous emission. 

In the discussed experiment (Robert et ai, 1992) a two-peaked hydrogen wave packet was prepared in an atomic 
interferometer in the internal metastable 2s i state, and arranged to enter a region of space with an electric field where 
the Stark mixing between the 2si and 2pi states led subsequently to spontaneous emission of a Lyman-a photon. 
Since the atom was moving, different parts of the wave packet were coupled to the field at different times and thus the 
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spontaneous decay was sequentially turned on. The nonsimultaneity of the coupling was not taken into account by 
Rzazcwski and Zakowicz (1992) which may be a source of the disagreement between the theory and the experiment. 

A nonsimultaneous coupling between the field and the atom will generally lead to a complicated, infinite set of 
coupled integro-differential equations that are quite difficult to solve even numerically. We have decided therefore to 
consider a problem which on one hand is simpler and on the other exhibits physical properties which are expected 
to occur if a nonsimultaneous coupling is involved. The results we discuss cannot be directly compared to the 
experimental data for two reasons. First of all, to simplify numerical computations we shall assume that the wave 
packet is a one-peaked and not a two-peaked Gaussian. Second, we model the nonsimultaneity by considering an 
ideal situation where the decay is turned on by a moving step function. The experiment was not that ideal and the 
sequential turn-on was much less abrupt. 

Nevertheless, in spite of these simplifications, the result we report is highly nontrivial. It shows that the "nonsimul- 
taneous spontaneous emmission" (NSSE) is drastically modified with respect to the SSE case. The radiation pattern 
exhibits a time-dependent anisotropy which depends on details of the turn-on. In the infinite velocity or infinite time 
limits it reduces to the ordinary Doppler- modified dipole pattern of Rzazewski et at. 



II. SPONTANEOUS EMISSION FROM A "SPACE-TIME SUPERPOSED" ATOMIC WAVE PACKET 

The standard analysis of SSE involves "equal-time and spacelike" superpositions of the atomic center-of-mass 
eigenstates. In what follows we shall consider superpositions of the center-of-mass eigenstates which additionally start 
to decay in a position-dependent way. One can say that the choice of the "i = 0" varies from point to point. We will 
concentrate on the simple case involving a single Gaussian whose decay is turned on by a propagating step function 
moving in the — z direction with a constant velocity v. Some preliminary results of such computations were briefly 
discussed in Czachor and Li (1996). More complicated examples, including the case of a "beaded" two Gaussian wave 
packet, will be discussed elsewhere. 

We begin with the nonrelativistic dipole (RWA) Hamiltonian 

H=^+H A + H F + H AF , (1) 

with 

Ha = 2^ tJ oO'3, 

k.fJ, 

H A f = - i Y^ n 9k^e %kr a+a klJl + h.c, (2) 

which are kinetic, internal, free field, and dipole interaction parts respectively. The dipole coupling strength is 
Qkti = £kp{e-p • £kfj,)/^, where e p is a unit vector of the dipole operator, and p is the dipole moment of the electronic 
transition; efc M , ^=1,2, are the photon polarization vectors, and £k = \/2Tthuj k /V is the electric field strength per 
photon in the quantization volume V. 

We assume the state of the atom-field system to be 

\i>s) = Jdra t (r)\r;+;0) 

+ I drP t {r,k,n)\r;-;k,n), (3) 

where 

/3t(r;fe ' M) = ( 2^ )l j d Ph& k ^Y (P ~ nkyr,h > ( 4 ) 

with p and r being the center-of-mass momentum and position, "+" ("— ") correspond to excited (ground) internal 
atomic state, and k, \i denote momentum and polarization of the emitted photon. 
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The solution of the above model for the case of SSE is given by (Rzazewski and Zakowicz, 1992) 

a t (p) = a (p)e- s « te (-*)-^+ s ») te « 



Ptfak,fj) = -«o(p)Sfc M " 



,-(7+«o)*e(<) 



S k t0(t) 



with 



so 



-Q(t)a (p)g ktJ ,- 



7 + s - s fe 
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(5) 
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2M 



- -HUJO + flUOk 



(6) 



where 7 determines the lifetime. The step function <d(t) is 1 for t > and otherwise. The 0(— t) in the exponents 
in (||) represents the fact that the evolution of the wave packet for t < is free and uncoupled to the vacuum 
field reservoir. To generate the NSSE solution we superpose the center-of-mass eigenstates with amplitudes (Q) 
but now instead of O(i), which is responsible for the simultaneous turn-on, we take a more general superposition 
with amplitudes containing 0(/(t,r)). The "switch-on function" f(t,r) depends on details of the experimental 
arrangement. Consider 
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(8) 
(9) 

(10) 



The state (|^) corresponds to the process where the emission is turned on by the step function 0(i + z/v) propagating 
in the —z direction with velocity v. (This should not be confused with a description in a moving frame!) Obiously 

/3 t ,oo{r,k,n) = p t (r,k,n) (11) 

as it should be. The t — > 00 limit of (0) will not differ from the analogous simultaneous case since the contributions 

from 9ml + z/(tv))^j for t — > 00 and v fixed are important only for the z — > 00 parts of the wave packet which are 

assumed to vanish. Nevertheless, this asymptotic behavior does not mean that the problem completely reduces to the 
simultaneous case. The transient modifications of the angular distribution shown in Fig. 2 are relatively long-lived 
and in principle can last arbitrarily long, depending on relative scales of r na turai, v and the size of the wavepacket. 

Denote t = r(t,z) = t + z/v. The probability density for NSSE of emitting a photon with momentum k and 
polarization [i is 



Vtkn = J dr\(3 t , v (r,k,fi)\ 2 



\9k^ 



2irh || Vjvs || 2 
exp 



dp J dp' z JdzQ{r) 



1 1 2 . 

! r, , . pz, — v y 



a o {p)a* (Px ,p y ,p' z )C(p)C* (p x ,p y ,p' z ), 



(12) 



with 
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£(p) = 



1 — exp 



o x k x 



Py k v i Pzk 
M M 



+ ^rr- - S k - h 



(13) 



where we have defined 8 k = tOk + ^/(2Mc 2 ) — o;q. Without loss of generality, we can put k y = 0, (which simply 
specifies the coordinates). In this work, we assume the wave packet is initially coherent and has a Gaussian distribution 
in momentum space, given by 



«o(p) = 



exp 



2 2 

pa 



In this case the integration over p y , p z and can be performed to obtain 



Vtkp, oc Jdp x J<IzQ(t) exp ^-2 7 r 
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where W(£) = exp(— £ 2 )erfc(— i£) is the exponentially scaled error function, and 
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The time-dependent width of the wave packet (due to spreading) is a 2 = a 2 — 2ihr/M. In the asymptotic limit of 
large t, only the second term in Eq. ([l6]) survives. Using the asymptotic expansion y/irr]W(r]) ~ i for r] — > oo, we can 
show that our general result for (|I^) reduces to the SSE case (Rzazewski and Zakowicz, 1992) if v — > oo. 
The total probability of spontaneous emission is now 



£>tfc P = Jdn k (l-\e p -e k \ 2 )P t (6), 

kf^i 

with Pt(#) the reduced probability distribution that is defined as 

P t (9) = Jdu k Q t (u k ), 
3 7 a 2 



(18) 
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and 9 is the polar angle between k and the z-axis. 

For the case of hydrogen 2pi — ► lsi transition, the parameters are 7 = (27r)50 MHz, Ao = 121.6 nm, the recoil 

velocity v re coii = hko/M = 3.25 (m/s) and energy w rec oil = hk 2 /(2M) — (27r)13.328 MHz. 

In Fig. [l], the results of the time-dependent (normalized) spectrum are shown for the case of v — 10(w reco ii) and 
v = l(w r ecoii). We have used a = Ao as normally assumed for a cold atomic wave packet in the recoil limit. With 
these parameters, the time a/v it takes to traverse the wave packet is much less than the lifetime of the metastable 
2si state, therefore neglecting its decay is a well justified approximation. We see that due to the time dependence 
introduced by the swept turning on process, a general broadening of the spectrum is seen, but in the limit t — > 00, 
it approaches the result known for the case of the SSE. This time-dependent broadening can be easily understood if 
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one considers a wave packet which is totally localized in a region of space where the step function has not yet arrived. 
The spectrum must then be completely flat (equal identically to 0) since the atom cannot radiate. However, after the 
exitation (or coupling) is completed the atom will radiate and the radiation pattern must be nonzero. Our results 
show that the pattern continuously deforms from a flat distribution to a Lorentzian. The dependence on the polar 
angle 6 exists as well but the results will not be explicitly shown here. Instead we show in Figs. |2[^ the total integrated 
signal (over the whole spectrum range) of the reduced (i.e. divided by the dipole pattern) probability distribution. 
This quantity should be constant in the case of a simultaneous spontaneous emission and is clearly modified by the 
non-simultaneity of "excitation" . The effect is transient and should asymptotically approach the dipole distribution 
for very large t. Note nevertheless that the range of time we consider at Fig. 2 is relatively large in atomic scales 
(— 50r na turai < t < +40r natura i in Fig. 2). For the wave packet we consider the effect should be observable for a few 
hundred of natural lifetimes and can be more long-lived for different wave packets and excitation velocities. 

The results suggest that when the time of transit of the "excitation" through the wave packet is of the order of 
the lifetime of the excited state, the properties of the spontaneous emission may be significantly modified. One may 
speculate that this kind of modification may be the source of the "fringes" observed by Robert et at (1992). 

To complete the analysis let us mention that we have taken the initial wave packet ao(r) to be completely coherent. 
A more general treatment would require treating both a*(r) and /3*(r) as statistical wave functions starting with a 
partially coherent a^{r) described by the spatial coherence function T(r, r') — (ao(r)ao(r')) . O ur result can be easily 
extended to the case of a partially coherent initial wave packet provided the coherence function T(r,r') is given. 
Depending on the particular form of T(r,r'), the properties of the spontaneous emission as given by the solution 
may be different. The richness of the structure we observe in such a simple case suggests that the approach 
developed here may allow for a determination of the partial coherence function T(r, r') by studying its NSSE. In a 
wider perspective, we can speculate that this relation is not only restricted to the case of spontaneous emissions. It 
should also be true for any other kind of non-simultaneous interactions, such as scattering of light from different parts 
of a wave packet at different times. In view of the recent success in achieving the Bose-Einstein condensation of alkali 
atoms (Anderson et al., 1995; Davies et al, 1995; Andrews et al. 1997; Mewes et al., 1997) we expect such a kind of 
non-simultaneous interaction scheme could provide a possible diagnostic scheme for detecting the partial coherence 
function, and revealing the off diagonal long range order of the condensate. 

III. FINAL REMARKS 

A closer look at the SSE solutions given in Rzazewski and Zakowicz (1992) shows that they very well agree with an 
intuitive picture of what is going on: We simply obtain a coherent superposition of the atom-photon states representing 
Doppler shifted emitted photons. The probability distribution of the atomic momenta manifests itself in an analogous 
distribution of the frequency shifts with respect to the dipole pattern. This explains the oscillatory behaviour of the 
spectral density of radiation which is completely unrelated to any photon interference. This is consistent with what 
we know about entanglement of quantum states and its role for a two-particle interference (Czachor, 1991). What is 
encourageing, however, is the fact that one could obtain an identical result without any center-of-mass calculations 
but simply by taking a superposition of suitably Doppler shifted Weisskopf-Wigner wave functions. 

This observation motivated the analysis presented in this Letter. We simply analogously superpose the Weisskopf- 
Wigner wave functions but in addition to the Doppler and recoil shifts we take into account the fact that parts of the 
wave packet which are "excited" at different times must have the initial time-of-decay "t = 0" different for different 
points of the wave packet. 

One may wonder whether on general grounds the approach is acceptable. For example, the moving step-function- 
excitation involves an interaction of the initial wave packet with some external fields. As such, it may lead to a 
decoherence of the atomic state. The decoherence-recoherence experiment reported by the MIT group (Chapman et 
al, 1995) shows however that an atomic wave packet may remain coherent after having interacted with a laser beam 
and even having spontaneously emitted a photon provided one afterwards purifies the state in a suitable correlation 
experiment. 

To close these remarks let us note that even the SSE is simultaneous in only one reference frame. And even though 
the nonsimultaneity experienced by a moving observer is not exactly of the kind we discuss in this Letter, it shows that 
superpositions of states corresponding to "different positions and different times of events" are physically unavoidable. 
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FIG. 1. Normalized spectrum (8 — 0) at three different times corresponding to the edge of the moving theta function at 
2 = 0.5Ao (dotted line), z — (solid line), and z — — 0.5Ao (dashed line). Dash-dotted lines show the unmodified Lorentzian 
natural lineshape. The step function moves along the z axis with velocity v in the negative direction (from right to left). 

FIG. 2. The reduced angular distribution for v = l(urecoii) with other parameters the same as in Fig. 1. as a function of time 
and 6. At t = the step function arrives at the center of the Gaussian. An analogous plot is completely flat in the SSE case. 

FIG. 3. The reduced angular distribution of radiation at the time corresponding to the edge of the moving theta function 
at z = — 0.5Ao for three different velocities (in units of u re coii)- The rest of parameters the same as in Fig. 1. For v = v TCCoa 
the atom tends to radiate in the direction of propagation of the "excitation". For v = O.lVrecoil the radiation is concentrated 
in the direction perpendicular to the excitation. With growing velocity of the excitation the distribution approaches the dipole 
pattern. 
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